Nuclear Multifragmentation in the Non-extensive Statistics - Canonical Formulation 
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We apply the canonical quantum statistical model of nuclear multifragmentation generalized in the 
framework of recently proposed Tsallis non-extensive thermostatistics for the description of nuclear 
multifragmentation process. The test calculation in the system with A — 197 nucleons show strong 
modification of the 'critical' behaviour associated with the nuclear liquid-gas phase transition for 
small deviations from the conventional Boltzmann-Gibbs statistical mechanics. 
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Most of the fragmenting systems are characterized by 
strongly off-equilibrium processes which cease due to dis- 
sipation. The theoretical description of the fragmen- 
tation process depends on whether the equilibrium has 
been reached before system starts fragmenting. If the 
equilibrium is attained, then the thermodynamic models 
using different statistical ensembles in a given fixed vol- 
ume (freeze-out volume) can be applied. The ingredients 
specific for the considered phenomenology enter through 
the definition of fragments sizes and (binding) energies, 
fragment internal excitation properties, system size, con- 
served quantities in this process etc. The observable 
characteristics of the fragmenting system are employed 
to fix certain features of, in general unknown, intermedi- 
ate equilibrium state. In nuclear physics, several models 
of this kind have been tried with unquestionable suc- 
cess in describing the transitional phenomenon in heavy 
ion collisions from the regime of particle evaporation at 
lower excitation energies to the explosion at about 5 - 
10 MeV/nucleon of the hot source accompanied by the 
copious production of the intermediate mass fragments 
jfl-S]. The situation when the fragment production has 
to be considered as an off-equilibrium process is described 
by various kinetic equations, mainly on the level of one- 
body distribution functions Q . Here the statistical equi- 
librium is not assumed but the kinetic models in turn are 
plagued by the unsurmountable conceptual difficulties in 
the calculation of asymptotic, observable features of the 
fragments. As an attempt to overcome at least some of 
these difficulties in both groups of models, in this work 
we extend the thermodynamic (canonical) model of the 
fragmentation in the framework of the recently proposed 
thermostatistics [^[ to include certain off-equilibrium cor- 
relations in the system. The Tsallis' generalized statis- 
tical mechanics (TGSM), which provides the basis for 
generating this new model, is based on an alternative 
definition for the equilibrium entropy of a system whose 
ith microscopic state has probability pi : 



tropy S q has the usual properties of positivity, equiprob- 
ability, concavity and irreversibility, and preserves the 
Legendre transformations structure of thermodynamics. 
In the limit : q — > 1, one obtains the usual Boltzmann- 
Gibbs formulation of the statistical mechanics. The main 
difference between the Boltzmann-Gibbs formulation and 
the TGSM lies in the nonadditivity of the entropy. In- 
deed, for two independent subsystems A, B, such that the 
probability of A + B is factorized into : pa+b — PaPb j 
the global entropy verifies : 

S q (A + B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B) (2) 

TGSM provides a natural framework for the thermody- 
namical formalism of the anomalous diffusion and ubiq- 
uity of Levy distributions || . Long-range correlations in 
the system, as appearing in the situation of the thermal- 
ization of a hot gas penetrating in a cold gas in the pres- 
ence of long range interactions, is typical for q > 1 In 
some cases, the entropic index q in TGSM can can be also 
related to the fluctuations of the temperature in the sys- 
tem ||. Variety of the off-equilibrium situations which 
can be accounted for within the TGSM make it useful 
as a basis for the generalization of the thermodynami- 
cal fragmentation models and , in particular in address- 
ing the problem of the influence of these non-extensivity 
correlations on the signatures of 'criticality' in finite sys- 
tems. In the context of nuclear multifragmentation, this 
is usually referred to as the signatures of liquid-gas phase 
transition in small systems. 

Our starting point is the canonical multifragmenta- 
tion model |9| with the recurrence equation method 
P,fio[ which makes the model solvable without the Monte 
Carlo technique and transparent to the physical assump- 
tions and generalizations. Canonical ensemble method 
in TGSM was discussed in pT| . The main ingredient of 
the TGSM generalization of the canonical fragmentation 
model p|] is the expression for fragment partition function 
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(i) 



Z q (s,t) = J2^ + ^(s,t)}- 1 ^ 



(3) 



and q (entropic index) defines a particular statistics. En- where q\ = q — 1, s and t are the fragment mass and 
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charge respectively, and the fragment partition probabil- 
ity equals : 



(4) 



where £j?(s,i) = p 2 /2m + U(s,t), f3 = 1/T. The inter- 
nal energy U includes the fragment binding energy, the 
excitation energy and the Coulomb interaction between 
fragments in the Wigner-Seitz approximation Q . Chang- 
ing summation into the integration in (||) one gets : 



z q ( s ,t) = (qr 1 -3/2)1^ (qr 1 ) 

X q- 3/2 [l + qi f3U( Sl t)]-^ +3/2 



(5) 



where g is the spin degeneracy factor, Vf is the free vol- 
ume and At = [(27r)/(mT)] 1 / 2 , where m is the mass of 
the fragment (s, i). In the limit q — > 1, one recovers the 
familiar expression : Zi = gVf\x~ 3 exp {—(3U). 

Given the partition function, the mean value of any 
quantity in TGSM is § : 



(0) 



For the average energy of fragment (s,t) one obtains : 



< e(s,t) > Q 



d fi-[z q (s,t)y 
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(7) 



In the dilute gas approximation |12| , a partition func- 
tion of the whole system can be written as follows : 



Q q (A,Z) = n 



[z q M] N ^ l) 

JVft («,*)! 



(8) 



where the sum runs over the ensemble Ha,z of differ- 
ent partitions of A and Z of the decaying system : 
{n} = {JV A (l,0),JV A (l,l),...,JV ft (i4,Z)} and N*(a,t) is 
the number of fragments of mass s and charge t in the 
partition n. In this approximation, the recurrence rela- 
tion technique fj] can be applied providing exact expres- 
sion for Qq(A, Z) : 



Q g (A,Z) 
1 
A 
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s,t<s;s-t<A~Z 



sZ q (s,t)Q q (A-s,Z-t) 



(9) 



These relations can now be conveniently used to calculate 
ensemble averaged characteristics. However, in order to 
ensure the proper normalization, it is better to work with 
generalized averages ]l3| : 



< O >g =< O >g I < 1 > Q 



(10) 



These normalized mean values exhibit all convenient 
properties of the original mean values (o). Moreover, 



when the normalized mean values ( |To| ) are used, the 
TGSM can be reformulated in terms of ordinary linear 
mean values calculated for the renormalized entropic in- 
dex : q* = 1 + (q — l)/q. In particular, the total average 
energy of the system becomes : 



S q = ^2< N(s,t) > q *az< e{s,t) > q , 
where < e(s,t) > g » is given in (0) and : 

< N(s, t) > qAZ = Z q (s, t) QM Q ~ { A Z z) t] 



(11) 



(12) 



Analogously, the heat capacity (— d£ q /dT \y) is : 
C v = P 2 {^2 E < A ( s t;s't) > q ,< s(s,t) > r x 

s .< s' ,t' 

x < e(s ,t) > g * +^2< N(s,t) > rA z x 

s,t 

x [<e 2 {s,t) > q , -<e{s,t) >%}} , (13) 

where : 

< A(st;st) > q =< N(s,t)N(s\t) > qAZ - 

- < N(s,t) > qA z< N(s',t) > qA z 



and 



< N(s,t)N(s ,t) > qAZ 
Z q (s,t)Z q (s ,t ) 

+ 6 ss ,5 tt ,Z q (s,t) 



Q q (A-s-s ,Z-t-t) 
Q q (A,Z) 
Q q (A-s,Z-t) 



Q q {A,z) 



(14) 



Fig. |l| shows the caloric curve for different values of 
the entropic index q > 1 in the above described canonical 
multifragmentation model. The calculations are done for 
the system with Z = 79 protons and N = 118 neutrons 
(A = 197). The free volume is : V f = 3A/ p = A/p f , 
where p = 0.168fm~ 3 . The excitation energy is : E* = 
£ g (T, pf ) — £g(T = 0, p ). Curves T{E*/A) for different q 
are very similar outside of the 'critical zone' of excitation 
energies : E*/A G (2.5 - 10)MeV. On the other hand, a 
strong sensitivity to even tiny changes of q can be seen 
inside of the 'critical zone'. 

This fragility of equilibrium (q — 1) 'critical behavior' 
to the small changes of the entropic index can be seen 
even better in Fig. || which shows the specific heat vs 
the temperature. Increasing of the entropic index q is 
associated with both a significant sharpening of the peak 
in Cy and an increase of the 'critical' temperature Tq. 
Since this increase of Tq is accompanied by only a small 
change of the total excitation energy (see Fig. Q), there- 
fore the kinetic part in the total energy increases with 
q. In other words, the multifragmentation in statistical 
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systems with q > 1 takes place in the hotter environment 
than in the limiting equilibrium case q = I. 
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FIG. 1. The 'caloric curve' for the system with A = 197 
nucleons is plotted for various values of the entropic index q 
: 1 (the solid line), 1.0005 (the dashed line)and 1.001 (the 
dotted line). 
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FIG. 2. The specific heat for the system with A = 197 
nucleons is plotted for various values of q : 1 (the solid line), 
1.0005 (the dashed line)and 1.001 (the dotted line). 

It is an open question whether the correlations for 
q =/= 1 change the nature of the equilibrium 'phase transi- 
tion'. Whereas the liquid-gas phase transition is char- 
acterized by properties of the largest cluster Q, the 
shattering phase transition in off-equilibrium systems is 



characterized by the multiplicity of fragments JTHj . Fig. || 
shows the average multiplicity dependence of the normal- 
ized second factorial cumulant moment of the multiplicity 
distribution : 72 = (< m(m — 1) > — < m > 2 )/ < m > 2 . 
72 is a measure of the fragment - fragment correlations 
and equals for the Poisson distribution. One can see 
strong build up of multiplicity fluctuations with increas- 
ing q in the 'critical region'. This enhancement of 72 is 
associated with the strong pic in Cy as seen in Fig. ||. 
For q = 1.001, the maximum of 72 is comparable with 
those found in the 2D and 3D percolation systems of 
comparable size p6|. 
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FIG. 3. The normalized second cumulant factorial moment 
of the fragment multiplicity distribution is plotted as a func- 
tion of the average multiplicity for various entropic indices q 
: 1 (the solid line), 1.0005 (the dashed line)and 1.001 (the 
dotted line). 

The fragment-size distributions dN/dA at T ~ Tc for 
different q values are shown in Fig. ^. One can see the 
significant evolution of dN/dA with increasing entropic 
index which, together with the evolution of multiplicity 
distributions (see Fig. ||), illustrate the change of mech- 
anism of the multifragmentation. At T = Tc, one hnds 
approximately power-like fragment-size distribution for 
q = 1, and the persistence of the heaviest residue for 
q > 1. The 'critical zone' for q > 1 is associated with 
the exponential fragment-size distribution. This resem- 
bles the critical binary fragmentation with the Gaussian 
dissipation ||l7j| which is an off-equilibrium process. For 
T/Tc > 1 and q > 1, the heavy residue explodes into the 
large number of light fragments and the fragment-size 
distribution remains exponential. 

In conclusion, we have developed the generalization of 
canonical multifragmentation model in TGSM . This new 
model provides an alternative way of taking into account 
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expected deviations from the thermodynamical equilib- 
rium due to nonextensive correlations in the multifrag- 
mentation process. We see the main advantage of the 
proposed approach in the correct description of of pro- 
duced fragments and in the preservation of mathematical 
structure of thermodynamics. Variability of different sig- 
nals of equilibrium phase transiton with respect to even 
small deviations from 5=1, demonstrates that the char- 
acterization of nuclear phase transition in terms of finite- 
size scaling analysis may be hazardous. On the other 
hand, signals of 'criticality' in (/-generalized canonical 
fragmentation model are stronger, what in turn demon- 
strates the robustness of this 'critical phenomenon'. It 
is worth mentioning that all considered characteristics 
change qualitatively with q in the 'critical zone' of exci- 
tation energies and one cannot exclude that order of the 
phase transition changes as well. 
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FIG. 4. The fragment-size (mass) distribution for the sys- 
tem with A = 197 nucleons (Z = 79 and N = 118) is plotted 
for entropic indices and different temperatures normalized by 
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The new flexible family of fragmentation models obey- 
ing q-statistics provides a powerful tool in analyzing ex- 
perimental data and in characterizing possible deviations 
from idealized equilibrium phase-transition picture in nu- 
clear multifragmentation in terms of the entropic index. 
This analysis is now in progress. 



4 



